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Abstract 

The possibility of repulsive Casimir forces between small metal spheres and a dielectric 

half-space is discussed. We treat a model in which the spheres have a dielectric function 
(N 

given by the Drude model, and the radius of the sphere is small compared to the correspond- 



ing plasma wavelength. The half-space is also described by the same model, but with a 
different plasma frequency. We find that in the retarded limit, the force is quasi-oscillatory 



This leads to the prediction of stable equilibrium points at which the sphere could levitate 
in the Earth's gravitational field. This seems to lead to the possibility of an experimental 



test of the model. The effects of finite temperature on the force are also studied, and found 
to be rather small at room temperature. However, thermally activated transitions between 
equilibrium points could be significant at room temperature. 
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I. INTRODUCTION 

The interaction energy between a particle with non-dispersive polarizability «o 
and a perfectly reflecting wall was found by Casimir and Polder |I[ to be 1 

VcP = "32^' (1) 

where z is the separation. One can assign a frequency spectrum <j[uj) to this potential 
and write |2| 

poo 

Vcp = i&hj duja(uj) ' (2) 

where 

(?{uj) = (2u 2 z 2 — 1) sin 2uz + 2ujz cos 2uz. (3) 

The integrand, ct(uj), is an oscillatory function whose amplitude increases with in- 
creasing frequency. Nonetheless, the integral can be performed using a convergence 
factor (e.g., insert a factor of e -7 ^ and then let 7 — > after integration). The result is 
the right hand side of Eq. (0). It is clear that massive cancellations have occurred, and 
that the area under an oscillation peak can be much greater in magnitude than the 
final result. This raises the possibility of tampering with this delicate cancellation, 
and dramatically altering the magnitude and sign of the force. 

In an earlier paper |2|, henceforth I, this possibility was investigated for small 
dielectric spheres near a perfectly reflecting wall. There it was found that the force 
can be considerably larger than that given by Eq. (Q), and furthermore can be either 
attractive or repulsive. The force was found to be quasi-oscillatory, with a period of 
the order of the plasma wavelength of the material in the sphere. This effect can be 
understood as a resonance phenomenon involving the vacuum modes and the plasma 
frequency of the sphere. The calculations in I were restricted to the case of a perfectly 
reflecting wall at zero temperature. In the present paper, we generalize these results 
to the cases when the wall has finite reflectivity and the temperature is non-zero. 

In Sect, m a force on a small sphere near a dielectric wall is computed at zero 
temperature. In Sect. 11111 finite temperature corrections are discussed. The results 
of the papers are summarized in Sect. II VI 



II. FORCE ON A SMALL SPHERE NEAR A DIELECTRIC WALL 

A small sphere of radius a is placed a distance z from a wall. Both the sphere 
and the wall are composed of uniform material characterized with dielectric function 



Lorentz-Heaviside units with c = U = 1 will be used here. 



e(u>). We will take the dielectric function to be that of the Drude model, 

e(u) = 1 - , Up , , , (4) 

where u p is the plasma frequency and 7 is the damping parameter. From now on, 
we take uo p to be the plasma frequency of the material in the sphere, and uo q to be 
that of the interface. We will use subscripts w and s to distinguish between other 
parameters describing the wall and the sphere respectively. At present, we assume 
that the temperature is zero. 

The mean Casimir force on the sphere can be written as a Fourier transform: 

F = i-|^FH, (5) 

where F(u) can be viewed as the contribution of a vacuum mode of frequency u. This 
contribution can be found from essentially classical considerations, as it is the normal 
component of the classical force on the sphere when a plane wave of frequency u is 
incident on the sphere-interface system. In I, it was calculated in an electric dipole 
approximation for the wave scattered by the sphere; the force is found as an integral of 
the Maxwell stress tensor over the sphere. In the present case, the extra complications 
come from the finite reflectivity of the interface. However, these can be handled using 
the dyadic Green's function techniques of Schwinger et al |3|. 
The result for F (or V) follows from 

FM = -Jv« r M(^, (6) 

where (E 2 ) is the (renormalized) expectation value of the square of the electric field 
at the sphere's location and a r (u) is the real part of the dynamic polarizability. This 
expression is equivalent to the familiar result that the interaction energy of an induced 
dipole with a static electric field E is 

V = -\aoEl (7) 

where «o is the static polarizability of the particle. It turns out that this expression 
can be applied to the dynamic case as well, e.g. by expressing (E 2 )^ in Eq. (0) as a 
transverse spatial Fourier transform 



[E\ = -Re {i J(dk ± ) j^l [cA + (2k 2 - u 2 ) r>] 



e -2nz 



where kj_ is the transverse wavevector, and k = |kj_|. The reflection coefficients due 
to two polarization states of the electric field vector, r and r', are given by: 



K — Ki 



K + K\ 
K€ w / 

K€ w -r K\ 



(9) 



r = ■ , (10) 



where k is defined by k 2 = k 2 — uj 2 , and K\ by k 2 = k 2 — uj 2 t w . 

A detailed discussion is given in the Appendix, but the result for the interaction 
potential (F = — VV) can be written as 

V = ReS.— f°° duja r f°° dkk—— [uj 2 r + (2k 2 - uj 2 ) r'] e~ 2KZ \ . (11) 



This is similar to Eq. (3.34) in Ref. [3j, which gives the interaction energy between 
a molecule and a dielectric plate. In our case, the molecule has been replaced by 
a dielectric sphere, and the polarizability of the sphere a has been replaced by its 
real part, a r . This is also a generalization of Eq. (45) in I. The treatment in I 
assumed no evanescent modes. However, the quantization of the electromagnetic 
field in the presence of a dissipative dielectric requires one to treat the frequency uj 
and wave number k as independent variables, effectively leading to evanescent wave 
contributions. A detailed justification for the use of the real part of the dynamic 
polarizability, a r , was given in I. 

The complex polarizability of the sphere is given by 

«H = «o^4^ (I 2 ) 

e s (uj) + 2 

where ao, the static polarizability, is given by «o — 4-7ra 3 . If e s (uj) is given by the 
Drude model, Eq. (0J), then the real part of the polarizability is 

uj 2 p -3uj 2 , N 

a r = a oj — (13) 

v (3u z — uj^Y + 9u Zr y; 

This function has four poles in the complex tu-plane, at ±Q ± «7 s /2, where 



^ = ^V ' l H-^s- ( 14 ) 

It will be convenient to deform the contour of integration in the cj-plane, and 
isolate the residues of these poles. However, we must first consider the location of 
other possible singularities of the integrand of Eq. (|llj). There are branch points at 
the values of uo for which k = and k± = 0. The former occur at uj = ±k. The 
latter are in the lower half u-plane. In the limit that 7^ <^C u q , they are located at 
approximately 

w=±v ^_ i ^L y . (15 , 

The dielectric function e w (uj) has a pole at u> — —ijw, but both of the reflection 
coefficients, r and r', are regular at this point. Finally, there is a possibility of poles 
in the reflection coefficients at points at which k = —K\ or en = —K\. However, it 
may be shown that no such points exist. The electric field Green's function should be 
defined by an integration contour which goes beneath the singularities for Re a; < 



FIG. 1: The contour of integration for Eq. (|11|) in the complex w-plane for fixed k is 
illustrated. Here X marks a branch point, and a dotted line is its associated branch cut. 
The poles of the function a r (ui) are marked with the + symbol. When the contour is rotated 
to the imaginary axis, the poles in the first and third quadrants are the only singularities 
encountered. 



and above them for Re a; > 0, that is, Feynman boundary conditions. However, this 
does not include the poles of a r (u), which is not part of the electric field Green's 
function. Thus the contour of integration is as illustrated in Fig. GJ 

We may now rotate the contour of integration to the imaginary cu-axis, so that 
V can be written as the sum of an integral over imaginary frequencies, Vj, and a 
contribution, V p , from the residues of the pole in a r : 



V = Vj + v p , 



with 



and 



V; 



i r°° r°° i 

— / d(a r (i() / dkk- [-( 2 r + (2k 2 + C 2 ) P] e" 
StH Jo Jo « 



2kz 



v„ 



p Re / dkk- 



48nQ J k(ui) 

where iO\ is the pole of a r in the first quadrant 



[w?f(wi) + [2k 2 - uo\) r'M] e - 2 ^ )z , 



U>! 



-^l2ui-9^ + -i ls , 



(16) 



(17) 



(18) 



(19) 



and r is denned as 



(20) 

and similarly for r'. The first term in Eq. (J2UJ) arises from the pole in the first 
quadrant, and the second term from that in the third quadrant, and we have used 
the fact that the Drude model dielectric function satisfies e(— u,j) = e(u, —7). 

By introducing polar coordinates u and ( ( = ucos9, k = usm9 ), and subse- 
quently taking t = cos(6 l ) , Eq. (JT7J) becomes 

(21) 
The numerical evaluation of the integral is done in the limit 7 S — > 0, 7^ — > (plasma 
model regime). In this case we can write the coefficients r and r' (in place of r and 
r') , as defined in Eq. (fTU|) . in terms of the new variables as 



u - Ju 2 + u 2 

v q 22 



u+ ^Ju 2 + ul 



u 2 t 2 + u 2 -ut 2 Ju 2 + u' 2 
r ' — 1 v 1 (23) 

n, 2-1-2 _|_ 1 *J, _|_ ni-f-2 //)l2 _|_ / .2 



y 



In Figures 121- IH the force on the sphere (F = —V (Vj + V p )) in units u^ a is given 
for three different ratios uo q /uo p . The V p term is typically dominant, and gives rise 
to a quasi-oscillatory contribution to the force, as seen in the figures. This quasi- 
oscillatory behavior is in contrast to the monotonic attractive force in the case of 
an atom in its ground state interacting with an interface |l|, [3J]. However, it is 
similar to the force between an atom in an excited state and an interface [6j], which is 
also quasi-oscillatory In the present case, the oscillations arise from a resonance at 
lu = uj\ ~ V3lu p /3, if 7 S <^C u p , as seen from Eq. (fI5j). The appearance of the real part 
of the polarizability is crucial for this feature of the result. The force is alternatively 
attractive and repulsive as a function of separation, and this leads to a sequence 
of stable equilibrium positions for the sphere. These are the zeros of the function 
F(z) with negative slope. This seems to lead to the possibility of an experimental 
test of the model in which the sphere could levitate in the Earth's gravitational field 
at positions where F cancels the sphere's weight, as discussed in I. Note that F is 
defined so that F > is repulsion, and F < is attraction. Also note that because 
the present paper adopts Lorentz-heaviside units, in contrast to the Gaussian units 
used in I, the scale in the figures differs by a factor of Air from Fig. 7 in I. 

Of particular interest is the effect of the finite reflectivity of the interface, which is 
described by its plasma frequency, u q . As expected, when u q ^> u p , the results agree 





- / Y 












— - perfectly reflecting wall 1 


1 




0.5 








V 















li 
It 
i 

i 






/,' 


V / 




\ I 
» /' 
V ' 




\\ 


s />' 






-0.5 




\ ji 
\\ I' 













2 4 6 



10 12 14 16 18 20 



FIG. 2: The force, in units of uj^ cto, is given for the case uj q = 5uj p . Repulsion corresponds 
to F > 0. Here the result is quite close to that for the perfectly reflecting wall (dashed 
line). 



with those given in I for a perfectly reflecting wall. As the ratio uj q /uj p decreases, the 
force oscillations shift their phase and eventually decrease in amplitude, as seen in 
the figures. In Figure El we see that if u q is larger than a few times u> p , the perfectly 
conducting result is a very good approximation (dashed line). Even when uo q < u p , as 
in Figure EJ one still finds the quasi-oscillatory behavior of the force. The amplitude 
of the force decreases, as expected, since the effect should disappear as u g — > 0. 

In the small distance limit, a -C z <^C ^ p l , the expression for the total force on the 
particle in the plasma model regime becomes 



l« + ^ 



tt'o 



3v/2 



UJ 2 p UJ q 



4vr 8 



2u 2 
p 



3uj 2 „z 4 ' 



(24) 



We see that the force diverges near the interface, and moreover its sign depends 
on the ratio u p /uj q , which seems to be an artifact of the assumption of a perfectly 
smooth interface. As discussed in Ref. |4J, dispersion alone is not sufficient to render 
mean squared electromagnetic fields finite in the limit that one approaches such an 
interface. 



A. Perfectly reflecting Wall 



In the limit e w — > oo, as seen from Q and (fTU|) . r 
case (I21J) can be analytically integrated over t yielding: 



-1 , and r' — > 1 . In this 
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FIG. 3: The force is given for the case cv q = 2u p . Here the finite reflectivity has caused a 
shift in the locations of the force maxima and minima, but little change of the magnitude 
of the maximum force. 



V; 



"A 



tt'o 



4tt 2 



X 



arctan 



2 V3 y/&% - 3^ 

V3(-2u + 



duu z 



Is 



and (|18jl. after analytical evaluation, becomes: 



arctan 



V„ 



a uJp 1 



VM2u+js 
VH 2 - 3 7. 



(2cj 2 2 2 + 2iwiz - l) e 



2iwiz 



-2«2 



(25) 



(26) 



p 4vr 24fi z 3 
It can be shown that these results agree with the ones in I, namely we find that 



the total force on the sphere in this limit becomes: 

F = J + P, 
where 



(27) 



J 



-i v > 



u„ 



-an 



16tt 2 z 4 



(se + u§ (4^ 3 e 3 + 6z*e + ^ + 3) c „ 2zg 

(3£ 2 + ^ 2 ) 2 -9£ 2 7 2 



(2^ 
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FIG. 4: The force is given for the case ui q = 0.5oj p . Besides a shift in the locations of the 
force maxima and minima, the finite reflectivity has caused a reduction in the magnitude 
of the maximum force. 
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V 



3"f 2 z 2 — 67^ — 6) sin 2Vtz+ 



1927rfi z A 
(l2 7s fiV - 7s 3 z 3 + 12fiV - 37 S V - 67^ - 6) cos2fiz] 



(29) 



As seen here, the plasma model is a good approximation as long as ^ s z <C 1. Oth- 
erwise, Eq. (|29|) can yield significant corrections for the more distant equilibrium 
positions, but not for the first several peaks. We expect this to be true for the case 
of finite conductivity of the wall as well, as long as 7 S <C u> p , and 7^ <C cu q . 



III. FINITE TEMPERATURE CORRECTIONS 

In the case of nonzero temperature, Eqs. (J 17)) and (|18|) have to be modified. We 
write Eq. (|17|) as a Fourier series instead of Fourier transform by the substitution 

47r 2 n 2 



C'-£ 



P 2 



2tt ~* p *-< ' 

n=0 



(30) 




FIG. 5: The force is given for the case uj q = lo p and uj p (3 = 5. The effect of the finite 
temperature is a slight increase in the magnitude of the force. 

The prime is a reminder to count the n = term with half weight, and j3 = 1/kT. 
Thus, in the limit 7 S ,7 UI ~ 0: 



V, 






Re E W^ f ^r [ ~ Cr + (2p + C) r ' ] e ~" (31) 

n=0 6 S,n + UJ p JO K 



where r and r', using Eqs. (JHJ), (|TU|) and (jlj) can be written as: 



K 



v^+ 



lc-'h 



r = 






(32) 
(33) 



We modify Eq. (fTHj) by inserting a factor f 1 + e/3 ^_ 1 j to account for the thermal 

energy. This factor reflects the fact that at zero temperature, each mode has an energy 
of |o>; at finite temperature, there is an additional thermal energy of 1/ [eP^ — l). 
The result is 



V„ 



a uj 



24vrfi 



^Re 



3 /3wi _ i 



+ - 1 f°° dkk-^ [w?r(wi) + {2k 2 - u\) r'(ui)] e' 2K{ ^ )z 

(34) 

The electric dipole approximation requires that the radius of the sphere be small 

compared to the dominant wavelengths involved, so our results require that a -C (3. 
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So long as the temperature stays within this range of validity, the thermal correction 
enhances the magnitude of the force, with no shift in the positions of the force maxima 
and minima as illustrated in FigureElfor case u p f3 = 5. In most of the cases of interest, 
uj p is in the ultraviolet range, so u p (3 ^> 1 at room temperature. The temperature is 
also limited by the melting point of the material of the wall. 

Although the finite-temperature corrections to the force are small, there is another 
thermal effect which can be more significant. This is the possibility of thermally 
activated transitions between the stable equilibrium points, which become important 
when the required energy is of order kT . We can estimate this effect by noting that 
for the case u q = 2u p , this energy is approximately 

/•7/w,, 

W= fdz^lAxlQ- 3 a uj* (35) 



Ji/uJp 

which can be expressed as 

W = 2200K (-2-V f^-V = 280K f-^-V f^-V . (36) 

\20nmJ VIOeW \10nm) VIOeW V ; 

Note that lOeV^ = l/20nm, so for uo p rj lOeV^ the upper limit on the size of a sphere 
that would fit inside of a particular potential minimum is about 20nm. In this case, 
the thermal activation effect would be small at room temperature. However, for 
a smaller sphere, e.g., a = lOnm, thermal activation would be noticeable at room 
temperature. 



IV. SUMMARY 

We have investigated a model for the amplification of vacuum fluctuations, in 
which a resonance near the plasma frequency significantly increases the magnitude of 
the force compared to the non- dispersive case. In particular, we have examined the 
effects of finite reflectivity of the wall and of finite temperature. Finite reflectivity 
can decrease the magnitude of the effect, but so long as u q is of the order of, or larger 
than u p , the perfectly reflecting results are qualitatively correct. Finite temperature 
can enhance the magnitude of the effect, but with typical plasma frequencies in the 
ultraviolet range, the correction to the force at room temperature is extremely small. 
However, thermal activation can be significant, especially for smaller spheres. 

The key prediction of this model is a force that is alternatively attractive and 
repulsive as a function of separation. The ratio of the peak force to that of gravity 
on a sphere near an interface with uj q > uo p is approximately given by 

Fmax „ 7 / w p \ 4 (lfim\ fig I cm 



<2yt 



Fy 



27 f U P Y fh^A ( l 9/cm Z \ e -5( 7 ,/leV)(*/lMm) j (37) 



where p is the mass density of the sphere. This suggests a possibility of an experimen- 
tal test of the model in which the sphere could levitate in the Earth's gravitational 

11 



field at equilibrium positions. If forces whose magnitudes are large compared to those 
give by Eq. {I} with cto ~ Ana 3 are found, the result can be interpreted as a form of 
amplification of vacuum fluctuations. The spheres need to be small (aoo p ^$ 1), so ex- 
periments on metal sphere would have to involve spheres with radii in approximately 
the lOnm to lOOnm range. In the case of gold spheres, for example, (u p = 9eV), the 
sphere radius would have to be less than 20nm and preferably less than about lOnm. 
One might also ask if a similar amplification is possible for the Casimir force 
between perfectly reflecting parallel plates. The spectrum of the Casimir force in this 
case was discussed by Hacyan, et al [7J, and one of the present authors [8J], and found 
to be quasi-oscillatory, as in the case of the Casimir-Polder potential. However, as 
was shown by Lifshitz J9( , the force between a pair of dielectric half-spaces is always 
attractive and no larger in magnitude than the Casimir force. Thus upsetting the 
cancellation seems to be more difficult for half-spaces, and suggests that the small 
sphere approximation may be crucial, to obtain a quasi-oscillatory force. For larger 
objects, there may be a cancellation of the contributions of different spatial regions. 
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APPENDIX A: DERIVATION FOR THE FORCE ON A SMALL SPHERE 
NEAR A WALL 

We can obtain Eq. (jUJ) by using Eq. (11) in I as a starting point for the force on 
the sphere 2 : 

F* = y>d 3 W + Ipj&E* + ?( P xB) i , (Al) 

where p is the dipole moment associated with the sphere, and E and B are the mean 
electric and magnetic field vectors at the sphere's location. We take p to be linearly 
related to E: p =aE. Using the Maxwell equation, E = V x B, we get for the last 
term in Eq. (jAlj) : 

~ (p x B) 1 = ~a [(V x B) x B] J = ~a [- (d l B k ) B k + (d k B l ) B k ] . (A2) 



2 



This expression is obtained by integrating the Maxwell stress tensor over a spherical surface just 
outside the sphere. 
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Thus Eq. (IA1|) becomes: 



1 



F l = -a [2E j d j E i + E j d i E j + 2B j &B i - 2Byd i B j \ 



(A3) 



Now, we replace the field products with their appropriate expectation values for i = z. 
The expectation values of the electric and magnetic fields can be expressed through 
the Green's dyadic as in J3j: 



and 



h 



(£,(r)£ fc (r')>=i>(r,r», 



(By[r) -Bfc(r')) = tji m e knv {did n lu) 2 )r mp (r, r', to). 



(A4) 



(A5) 



Quantities such as (Ej(r) £^(r')) must be real, so we need to take a real part, which 

will only be done explicitly in the final expressions. Some components of T are 
(here k x is chosen to point along the +x axis): 



T xx (r,r',uj) 



dk. 



(27T) 5 



p ikx(r-r')j 



X 



r w (r, r» 



r M (r,r» 



i 19 19 / e -K|*-*l +r .' e -*(*+*) 

— 5(z - z') + -- — — 

e V ; edze'dz'\ 2k 

1 , , v /„-K|z-z'| _|_ r „-K(2 + z') 



dk 



1 



-a; 



2k 



,ik x (r-r')j 



X 



| ^2 /g-re|z— z'\ I r > e —K,(z+z') 



-JU-:') + - 



T xz (r,r',uj) = dk ± —— I e 



1 



3 ikj.(r-r')j 



2 fi- 



fe d ( er K \ z ~ z '\ 4- r * e -K(z+z') 



(27T) 2 

1 

{2*r 



(A6) 

(A7) 



ee' dz V 2k 

ee' dz' 



(A8) 
(A9) 



r, ( .(r. r'.^) = / r/k ± ^^(r-r'U (-») -^ ( — — - -^(AlO) 



2k 






Hence, for the first term in Eq. (jA3|) we have: 



i (E j dyE z ) = d x ,T xz + d y ,T yz + d z ,T zz . (All) 

Next we drop the (^-function terms, and take the coincidence limit, r' = r and z' = z, 
after performing the dy differentiation. All derivatives in y' are zero, since k^ points 
along the +x axis, so the second term above is zero. Using Eqs. (jA9|) and (jA8|) . we 
find the remaining terms (with e = e' = 1 for the vacuum region): 
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Q*V = -\fdk ± j^k 2 (r'e~ 2KZ ± 1) , (A12) 

d z ,Y** = -ljdk ± j±y 2 k 2 (r'e- 2 ^ T l) . (A13) 

Here the sign of the last term is determined by whether z approaches z' from above 
(upper sign) or from below (lower sign). We will argue later that these terms with 
ambiguous sign do not contribute to the final result. Now Eq. (jAll)) becomes: 

i (E j d f E z ) = - [dk ± -\,k 2 r'e- 2KZ . (A14) 

J {2n) z 

For the second term in Eq. ()A3)) we have: 



i {E 3 d z E>) = d z , (T xx + T yy + V zz ) . (A15) 

Using Eqs. (jMJ), fATjl . and (jMj) , we can write: 



d z ,T xx = ~\Jdk L -±-y 2 K 2 (r'e- 2KZ ± l) , 


(A16) 


d*r m = ~ljdk ± ^u 2 (re- 2 - T 1) , 


(A17) 


d*T„ = ~\fdk L j±y 2 k 2 (r'e- 2KZ T l) , 


(A18) 


so that Eq. (JA15J) becomes: 




i (E^BP) = l fdk ± l [(2k 2 u 2 ) r'e- 2KZ + co 2 re' 2KZ T 2cu 2 ] . 


(A19) 



For the third term in Eq. ()A3j) . we can write: 

'Bj&'B*} = (B x d x ,B z ) + (B y d y ,B z ) + (B z d z ,B z ) . (A20) 

The nonzero terms, using Eqs. (jA5|) and (jA7|) . become: 

i (B x d x ,B z ) = ~ljdk ± -±-k 2 (re~ 2KZ ± l) , (A21) 

i (B z d z ,B z ) = -\Jdk ± -^k 2 (re- 2 - T l) , (A22) 
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So, Eq. IjCTjl becomes: 



i(B^B') = ~ljdk ± j±y 2 k* 2re- 2 ^. (A23) 



To find the fourth term in Eq. (jA3|) . note that 

dk\ 



x 



„— k|z— z'| I „„— k(z+z') \ /„— re|z— z'l _i_ „/„— k(z+z') 

(fc 2 + V z V 2 ,)(- V^ Ww*' 



2k J V 2k 

From here we get: 



(A24) 



i (b^B*) = id* {B 3 B 3 ) = ~\Jdk± j^y 2 { [(2A; 2 -uo 2 )r + uo 2 r'] e~ 2KZ T 2lu 2 } . 

(A25) 
Combining Eqs. (lATIj) . (JA"19|) . (|A"23J) . and (lA"25J) . as in Eq. (|A"3J) . we have: 



2 (E j d f E z ) + (E^'EA + 2 (Bj&'B*) - 2 (B j d z 'B j 

Re (^/Wl jijj I 3 K 2 ^ - ^ 2 ) r' + ^ 2 r] e" 2 " ± 2^ 2 }) . (A26) 

The last term above would give an infinite force and cannot be present. If we average 
over z > z' and z < z', it would average to zero. In the limit that r = r' = 0, the 
force must vanish, so we can drop the last term. Then, we let a — ► a r , so that 
Eq. (|A3|) becomes: 



F» = Re iljdk ± Jr^M [{^k 2 - uj 2 ) r> + u 2 r] e~ 2 A (A27) 

This is equivalent to Eq. (0), with (E 2 )^ as defined in Eq. (jHJ). 
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